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Abstract 

We obtain a new solution of the star-triangle relation with positive Boltzmann weights which 
contains as special cases all continuous and discrete spin solutions of this relation, that were 
previously known. This new master solution defines an exactly solvable 2D lattice model of sta- 
tistical mechanics, which involves continuous spin variables, living on a circle, and contains two 
temperature-like parameters. If one of the these parameters approaches a root of unity (cor- 
responds to zero temperature), the spin variables freezes into discrete positions, equidistantly 
spaced on the circle. An absolute orientation of these positions on the circle slowly changes be- 
tween lattice sites by overall rotations. Allowed configurations of these rotations are described 
by classical discrete integrable equations, closely related to the famous (34-equations by Adler 
Bobenko and Suris. Fluctuations between degenerate ground states in the vicinity of zero tem- 
perature are described by a rather general integrable lattice model with discrete spin variables. 
In some simple special cases the latter reduces to the Kashiwara-Miwa and chiral Potts models. 
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1 Introduction 



There are only a few exactly solvable models in statistical mechanincs where the Yang-Baxter equa- 
tion takes its distinguished "star-triangular" form. The most notable discrete-spin models in this 
class are the Kashwara-Miwa [1] and chiral Potts [2Hlj models (both of them also contain the Ising 
model [5j and Fateev-Zamolodchikov ^Tv-model [6] as particular cases) see [7] for a review. There are 
also important continuous spin models, including Zamolodchikov's "fishing-net" model 0, which de- 
scribes certain planar Feynman diagrams in quantum field theory, and the Faddeev-Volkov model [9], 
connected with quantization [lOj of discrete conformal transformations [11^112]. All above models 
are also distiguished by the positivity of the Boltzmann weights — the property that is required for 
many applications, but rarely realized for generic solutions of the Yang-Baxter equation. 

In this paper we present new solutions of the star-triangle relation which possess the positivity 
property. Most importantly we present a master solution, which contains as special cases all the 
previously known continuous and discrete spin solutions mentioned abov^l], and also leads to new 
ones. Our master solution involves continouos real-valued spins, varying in the range < x,y < 2tt, 
and a reflection- symmetric Boltzmann weight, which is unchanged upon interchanging the spins x, y, 

<P(x — y — \a) ^{x + y — \a) 

and depends on an additive spectral parameter a (it enters additively into the Yang-Baxter equation 
p.5p below). Here is a spin-independent normalization factor and ^{s) is the elliptic F-function 

IMSll, 

Ep-ins I 
^T7 T > ■ (1-2) 

The latter depends on two fixed parameters p and q (elliptic nomes), 

p = e'''^ , q = e'''^, r] = -\tt {t + a) , Imr > 0, Imcr > . (1.3) 

This function obeys a simple functional equation <I>(s) *!?(— s) = 1, which ensures the reflection 
symmetry of the weight (|1.1|) . Define also a single-spin weight 

e^/4 

S{x) = _i?i(x|r)i?i(x|a) , (1.4) 
47r 

where 'djlz \ r), j = 1, 2, 3, 4, are the standard Jacobi theta-functions [16] with the periods vr and ttt. 
We state that the above weights satisfy the star-triangle relation of the form, 

dxo cS(xo) Wn-ai{xi,xo) Wo-i+asla^a, a^o) y^v-asi^s^^o) 

(1.5) 

— cti— as 

(xi,X3) Wag (Xi, 3:2), 



where the crossing parameter rj is defined in (jl.3p and IZ is some explicitly known scalar factor (see 
()2.34p below) which depends on the spectral variables ai and 03, but is independent of the spins 

Xi,X2,X3. 

^To be more precise, it only contains the solutions, which have a single one-dimensional spin at each lattice site. 
For this reason, it cannot contain the D > 2 fishing-net model which has multi-dimensional spins. 

■^Our definition (|1.2|l differs slightly from the standard definition of the elliptic F-function, see (|2.27|l below. The 
sum in (|1.2|) is taken over all positive and negative integers, excluding zero; it converges in the strip | Im s \ < rj, while 
the product formula (|2.24p applies to the whole complex plane of the argument s. 
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The weights and (|1.4p are real and positive in (at least) two main physical regimes 

(i) p = p*, q = q* or (ii) p = q* , (1.6) 

with real spins and a real spectral parameter in the range < a < rj. Note, also that the weights 
are unchanged upon negating the spins yVaix,y) = WQ,(±x,y) = Waix, ±y), S{x) = 5(— x), and 
are periodic in each spin argument 

Wa{x,y) = Wa{x + 2TT,y)=Waix,y + 2TT), 5(x) = 5(x + 27r), (1.7) 

therefore one can regard the spins as angle variables on a circle. This also means that the integral in 
p.Sp is a closed contour integral, where the integration contour can be deformed into the complex 
plane, if necessary. 

As is well known [T7] every solution of the star-triangle relation can be used to define exactly 
solvable edge interaction models on various two-dimensional lattices. For purposes of this intro- 
duction it is enough to consider an homogeneous square lattice. In this case the partition function 
reads 

2= I ■■■ I Wy^a{xi,Xj) ^Wr^-a{xk,Xi) ^S{Xm)dXm (1-8) 

where the first product is taken over all horizontal edges (ij), the second over all vertical edges 
{kl) and the third over all internal sites m of the lattice. We will implicitly assume fixed boundary 
conditions. In the limit of a large lattice the partition function can be calculated with the inversion 
relation method |18H20j : the result is given in (j2.42p . 

The elliptic F-function (jl.2p first appeared implicitly in Baxter's pioneering paper [21] on the 
8-vertex (it enters the exact expression for the partition function) and then was developed system- 
atically in [13H15j . Our key observation is that as a mathematical identity the star-triangle relation 
(jl.Sp reduces to the most general form of Spiridonov's celebrated elliptic beta integral [15]. Actually, 
it is quite remarkable that this fundamental integral identity, which lies at the basis of the theory of 
elliptic hypergeometric functions [22], is nothing but a Yang-Baxter (star-triangle) relation, defining 
a perfectly physical integrable lattice model of statistical mechanics. 

Let us now explain how this continuous-spin model could turn into a model with discrete spin 
variables, for instance, into the chiral Potts model. The key to the answer is the low-temperature 
limit. Note that the weights (jl.ip and (jl.4p symmetrically depend on two temperature-like parameter 
p and q. For generic values of these parameters, < |p|, |q| < 1, the formula (|l.ip defines a smooth 
function of spins with two dull bell-shaped maxima near x = y and x = 2tt — y. However, when one 
of the parameters p, q approaches the unit circle, the maxima become very sharp. Moreover, the 
function (jl.ip starts to exhibit additional sharp (5-function type) maxima, so that the spin variables 
become locked to a discrete set of energetically favourable position^. To get a better idea of how 
this happens, consider the limit when q approaches a root of unity 

p = e^^^ = fixed, q = e-''^e™/^, T ^ 0, N>1, (1.9) 

and p is some suitable chosen numerical coefficien10. The partition function (jl.8p develops a typical 
low-temperature asymptotics 

Z=j---j exp(-^ + 0(l) + 0(T)) T^^O, (1.10) 

0<Xm<'ilT ^ 

^Cf. a similar phenomenon for q — !> — 1 limit in ref. |23|. 

''Note that this limit does not belong to either of the main physical regimes (|1.6p . however, with a suitable choice 
of r and a it can be mapped to another physical regime for at least two terms of the of the low-temperature expansion 
discussed below. 
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where X = {xi,X2, ■ ■ ■ ,xm} and M is the number of internal (non-boundary) sites of the lattice. 
Explicit expressions for the energy functional £{^X^ are given in Section [3.31 and Appendix [Al It is 
a bounded from below function of the spin variables and that for > 1 it has stronger periodicity 
properties in each spin variable 

£[xi,X2, ...,Xm+'^ .,Xm) =£{xi,X2, . . . , Xm , ■ ■ ■,Xm), m = 1, . . . ,M , (1.11) 

than could be expected from (jl.7p . To obtain the leading asymptotics of the partition function (jl.lOp 
at T — )• 0, one has to minimize £ (A) , 

log^ = ^—-^+logZo + 0{T), (1.12) 

where 

d£{X) 



= 0, m = l,2,...,M, (1.13) 

and A^*^^) = {Ci; ^2, • • • , Cm} denotes an stationary point of the functional £{X), corresponding 
to the ground state of the system. Remembering that the spin variables run over the full circle 
< Xm < 27r, one immediately concludes that this ground state is A*^- fold degenerate. Indeed, 
thanks to (|l.lip . the equilibrium positions can be independently shifted, 

U^(m + 2TTnm/N, G Z^v = {0, 1, . . . , A - 1}, (1.14) 

without affecting the validity of the conditions (|1.13p . Thus at zero temperature, T = 0, the system 
"freezes" into one of these degenerate ground state configurations. 

The next natural question is to understand what is happening in the vicinity of the zero tem- 
perature, when the above degeneracy is lifted. A simple analysis shows that the only fluctuations, 
that are allowed in the next-to-leading order of the low-temperature expansion, are the discrete flips 
(I1.14P between different ground state conflgurations, whereas the values of ^m (mod 2tt/N) arising 
from the minimisation of £{X) remain frozen. This means that the quantity Zq in constant term 
of the expansion (I1.12P can be understood as the partition function of a certain edge interaction 
model with discrete spins variable {ni,n2, . . . jUm}, each taking A different values Uk E ^at. Its 
Boltzmann weights can be found by expanding (jl.Sp in the limit (|1.9p and isolating all constant 
terms. The results are given in Section [H Naturally, these Boltzmann weights explicitly depend on 
the original additive spectral variables a and r] — a, and the remaining temperature-like parameter p 
(which is still at our disposal). Moreover, they also depend on the variables {^mji which solves the 
minimization equation (|1.13p . and this is what makes the problem really complicated. In particular, 
the variables {Cm} could depend on the original spectral variables a and rj — a in a very complicated 
way, therefore, in general, the simple addition law for the spectral parameters in the star-triangle 
relation will be lost. Nevertheless, we definitely know that the emerging discrete spin model must 
be integrable! Indeed the original master model is integrable for any value of parameter q, therefore 
this integrability should manifests itself in every order of the expansion (I1.12P . 

It is clear, of course, that to analytically describe the new discrete spin model one needs to 
better understand the minimisation equations (jl.lSp . which according to our expectations must be 
integrable as well. The energy functional £{X) is a sum of two-spin edge energies therefore the 
variational equation (jl.l3p for any particular spin Xm will also involve spins on all neighbouring 
sites. For the square lattice each site has four (nearest) neighbours, therefore each equation (I1.13P 
will contain five variables. For the homogeneous model (jl.Sp on the square lattice they have the 
same form for any internal spin C, 
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Figure 1: The planar graph (shown by open circles and bold edges) and its medial graph ^ (shown by 

thin edges and alternatively shaded faces). 



where t,u^^ditt^^r are the spins immediately above, below, to the left and to the right of ^. The 
function ^'j(x,y) reads 

T?,(f (x - y + ia) I Nt) ^i(^(x + y + \a)\ Nt) 

^,(x,y) = — --^ j = 1,2,3,4- (1-16) 

(x -y-\a)\ Nt) (x + y-\a)\ Nt) 

Our next important observation is that these non-linear difference equations are indeed integrable, 
as expected. They appeared previously in remarkable papers by Adler, Bobenko and Suris |24tl25j 
devoted to the classification of classical integrable equations on planar quadrilateral graphs. More 
specifically, the minimization equation (|1.13|) arising here is closely related with the so-called Q4 
equation, which is the most complicated equation located at the top of the Adler-Bobenko-Suris 
classification. 

The organization of the paper is as follows. In Section 2 we briefly review some basic facts 
from the theory of integrable lattice models and present the master solution (jl.ip . (jl.4p . The low- 
temperature expansion is considered in Section 3. In Section 4 we present a new discrete spin solution 
of the star-triangle relation and show that in some simple special cases this solution reduces to those 
of the Kashiwara-Miwa [Ij and chiral Potts [2H3] models (the latter were the most complicated 
solvable discrete spin models hitherto known). Here we present only final results; the details of 
calculations will be published separately. 



2 The star-triangle relation 

2.1 Edge interaction models 

To facilitate further discussions let us briefly review some basic facts from the theory of integrable 
lattice models. A general solvable edge-interaction model on a planar graph can be defined in the 
following way [26^127]. Consider a planar graph of the type shown in FigUl where its sites (or 
vertices) are drawn by open circles and the edges by bold lines. The same figure also contains 
another graph shown by thin lines, which is the medial graph for ^. The faces of =Sf are shaded 
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alternatively; the sites of ^ are placed on the unshaded faces. We assume that for each line of one 
can assign a direction, so that all the lines head generally from the bottom of the graph to the top. 
They can go locally downwards, but there can be no closed directed paths in J£ . This means that 
one can always distort without changing its topology, so that the lines always head upward^. 

Now we define a statistical mechanical model on ^ . With each line I oi ^ associate its own 
"rapidity" variable "pi. At each site i oi'^ place a spin Sj, taking some set of (continuous or discrete) 
values. Two spins interact only if they are connected by an edge. This means that each edge is 
assigned with a Boltzmann weight which depends on spins at the ends of the edge. The Boltzmann 
weights usually depend on some global parameters of the model which are same for all edges (for 
instance, temperature-like variables). Here we assume that they also depend on local parameters, 
namely, on the two rapidity variables associated with the edge. 



> 



a O- 



■O h 



V 9 

M^P9(a;,2/) 



Figure 2: Edges of the first (left) and second types and their Boltzmann weights. 

The edges of ^ are either of the first type in Fig. [2l or the second. Let a, h be the spins at the 
end sites of the edge and p, q the rapidities of the associated line^EI, arranged as in Fig. [2l Then if 
the edge is of the first type, the spins a, h interact with the Boltzmann weight function VFpg(a, 6). If 
the edge is of the second type, they interact with the weight VFpg(a, 6). In general, there may also be 
a single-spin self-interaction with a rapidity-independent weight S{a) for each spin a. The partition 
function is defined as 

^ = E n^P5(«*'«i) V^v<ii^^^'^^ V^S{sm) (2.17) 

{s} (ij) (kl) 

where the first product is over all edges {ij) of the first type, the second is over all edges (kl) of the 
second type, and the third one is over all sites m. The sum is taken over the values of spins on the 
internal sites, while the boundary spins are kept fixed (for continuous spins the sum is replaced by 
an integral). The integrability requires the weights to satisfy the two star-triangle relation [7] 

^S{d)Wpg{d,c)Wpr(b,d)Wgria,d) = RpgrWpg{b,a)Wpria,c)Wgrib,c) 
d 

(2.18) 

Y,S{d)Wpg{c,d)Wpr{d,b)Wgrid,a) = RpgrWpg{a,b)Wpric,a)Wgric,b) 
d 

''This assumption puts some restrictions on the topology of the planar graph but still allows enough generality 
for our considerations here. 

''To avoid confusions note that the rapidity variables p and q are not related to the elliptic nomes p and q in H1.3[l . 
denoted by upright symbols. 
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where Rpqr is some factor independent of the spins a, b, c. For ah known solutions it can be written 
in the form [71|28], 

Rpqr — fqrfpq/ fpr-i (2-19) 

where fpq is a scalar function of the two rapidity variables p and q. Moreover, the weights can always 
be normalized so that they satisfy the two inversion relations, 



^Wpq{a,c)S{c)Wqp{c,h) = S{a) ^fpqfqpdab, 
c 

Wpq{a,b)Wqp{a,b) = 1, Va,6. 



(2.20) 



The first star-triangle relation in ()2.18p equates partition functions of the "star" and "triangle" 
graphs shown in Fig. [3l where the external spins a, b and c are being fixed. Similarly, the second 
relation in (I2.18P is related to a mirror image of Fig. [3l Note that the two relations coincide for 




— b 




Figure 3: A pictorial representation of the first star-triangle relation in (|2.18p . 

reflection-symmetric models, 

Wpqia, b) = Wpq{b, a), Wpq{a, b) = Wpq{b, a), (2.21) 

where the weights are unchanged by interchanging the spins a and b. 

The partition function (j2.17p possesses remarkable invariance properties |261l27j. It remains 
unchanged by continuously deforming the lines of ^ with their boundary positions kept fixed, as 
long as the graph ^ remains directed. It is easy to see that all such transformations reduce to a 
combination of the moves, corresponding to the star-triangle (j2.18p and inversion relations (12.20p . 
In general the partition function acquire simple fpq factors under these moves, however with an 
appropriate normalization of the Boltzmann weights the invariance is strict (all extra factors could 
be eliminated, see (I2.36P ). Given that the graphs ^ and ^ can undergo rather drastic changes, 
the above invariance, called "Z -invariance'' in [26], is rather non-trivial. It provides an ultimate 
formulation of the integrability statement. For instance, the commutativity of transfer matrices in 
solvable models is a particular case of this invariance. 

2.2 The master solution of the star-triangle relation 

The star-triangle relation (jl.5p is a particular case of (j2.18p . where the Boltzmann weights are 
reflection-symmetric and possess the difference property, which means that they depend only a 
difference of two rapidities p and g, 

Wpq{aM ^ Wp_g(a,6), Wpq{a,b) ^ W^-p+q{a,b). (2.22) 



7 



With this correspondence the spectral variables ai and 0:3 in p.5|) are connected to the rapidity 
variables in (j2.18p as 

ai= q-r, 03 = p - q . (2.23) 

Note also that the fact that the weights for the edges of two types in Fig. [2] are obtained from each 
other by a simple substitution a — t- r/ — a of the difference variable a = p — q is called the crossing 
symmetry. By this reason the parameter 77 in (j2.22p is usually called the "crossing parameter" . 
The standard elliptic F-function is defined as (see, e.g., Eq.(2.18) of [22]). 

, , 00 -1 2n+2 „2m+2 

n,m=0 ^ ^ 

where 

p = e''^^ , q = e""" , (2.25) 
and 0<|p|,|q| <1. Define the crossing parameter r/ by 

e"^ = pq , T] = -i7r(r + a). (2.26) 

For our purposes it is convenient to use the function 

Hx) F(e-'(^-"');p2,q2) = exp | V ^ , (2.27) 

[^on(p"-p-")(q"-q-") J 

which satisfies a simple "reflection" equation ^(x)^{—x) = 1. 

Our master solution of the star-triangle equation (jl.Sp (already quoted in the introduction) reads 

, , . X 1 ^(x — y + \a) ^(x + y + \a) 

Wa{x,y) = K(a)-i-) '—^ ) J . , (2.28) 

— y — \a) + y — \a) 

e'?/4 

S(x) = Mx\t)Mx\(t) ■ (2.29) 

where '&j{z \ r), j = 1, 2, 3, 4, are the standard Jacobi theta-functions [16] with the periods vr and vrr. 
The weights are symmetric, 

Wa{x,y) = Wa{y,x) = Wai±x,y) = Waix,±y), Six)=S{-x), (2.30) 

periodic in their spin arguments, 

Wa{x,y) = Wc.{x + 2Tr,y)=Wa{x,y + 2Tr), S{x) = S{x + lir), (2.31) 

and satisfy the recurrence relations 

Wa{x - TTa,y) '&i{^{x - y + \a)\T) -i?4(i(x + y + ia) | r) 

= , and the same with t a . 

Wa{x + Tra,y) i^^i^ix - y - \a)\T) i?4(^(a; + y - ia) | r) 

(2.32) 

For real-valued spins the weights are real and positive in the two main physical regimes 

(i) p = p*, q = q* or (ii) p = q* , (2.33) 
provided the spectral parameter is real and kept in the range < a < rj. 
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As a mathematical identity the star-triangle relation (|1.5|) contains seven continuous parameters: 
two nomes p and q; two spectral variables ai and a^; and three spins xi,X2,X3. One can show that 
this relation can be reduced to the most general form of Spiridonov's celebrated elliptic beta integral. 
Then using Eq.(3.1) of [22j one obtains the expression for the factor TZ in p.Sp . 

■^^fpll^, ;(„) = $(i,_2i„)^*!) , (2.34) 

/(ai + aa) k{7] - a) 

It is worth noting that if the function k(q) satisfies the equations 

^'^^ = «>(ir?-2ia), K(a)K(-a) = 1 , (2.35) 

K{a) 

then one trivially obtains 

f{a) = l, n=l. (2.36) 
With this normalization the inversion relations (|2.20p for the weights (j2.28p and (j2.29p simplify to 

2tt y 



dzS{z)Wr^-a{x,z)Wr^+a{z,y) = -^^^^^ {5{x - y) + 5{x + y)) 



(2.37) 



Wa{x,y)W^a{x,y) = 1, 
where 5{x) denotes the periodic 5-function 

^ oo 

K^) = ^ E e'"' • (2.38) 



27r 

n=— oo 



Moreover, for the same normalization the weights have the following special values 

Wa{x,y) =1, W^-a{x,y) = it^AK^ - v) + + y)) (2.39) 

0=0 a-!>+0 26 (x) 

Note that all the above relations are consistent with the symmetries of the weights ()2.30p . 
The general definition of the partition function (j2.17p specialized to this case leads to 

Z= j WWe^A^'^^Xj) J{S{Xra) dXm (2.40) 

(ij) 

where the edge variables are given by 

{p — q, for a first type edge 

r] — p + q, for a second type edge 



(2.41) 



where p and q are the pair of rapidities, associated the edge (ij), as shown in Fig. [2j The integral 
in (j2.40p is over all internal spins. The boundary spins are kept fixed (e.g., all boundary spins can 
be set to zero). 

The equations (|2.35p are precisely the functional equation of the inversion relation method |18fl2U] 
adapted for edge interaction models [29] . The quantity K{a) is interpreted there as the partition 
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function per edge in the limit of an infinitely large lattice. A suitable solution of (j2.35|) with 
appropriate analytic properties is given by 



K{a) = exp < 



If this function is used for the normalization of the Boltzmann weights (|2.28p in ()2.40p then the 
partition function per site 

lim Z^/*^ = 1, (2.43) 

M^oo 

is equal to one in the thermodynamic limit0. 

3 Low-temperature expansion 

3.1 Asymptotics of the weights 

Consider the limit of star-triangle equation (jl.Sp when one of the temperature-like parameters tends 
to a root of unity, 

p = e*^^ = fixed, q = e'^e'^/^ , £ ^ 0, N>1 . (3.44) 
Note that the crossing parameter (jl.3p then becomes 



% = -i(^+-r + ^)^^ = -i(^+.r) . (3.45) 

As explained in the Introduction the Boltzmann weights develop a singular asymptotics (ll.lOp where 
the leading term is unchanged upon the shifts Xi ^ Xi + 2tt/N . Therefore it is convenient to set 

From now on we will use new variables 

^' + ^ ' ^ = , : ' V = , 'r = . (3.47) 



1 + Nt 1 + NtV 2 J 1 + Nt' ' 1 + Nt 

instead of a, r and tj. With the new variables the asymptotics of (j2.28p and (j2.29p can be written 
a^ 

logW„(xi,X2) = -^£{e\^i,^2) + logWe{(l)i,(l)2,ni,n2)+0{e) , (3.48) 
log5(xi) = --C(0i)-^loge + log5((/.„n,) + O(e) (3.49) 



where 



WO = i(M' (3.50) 



'^Given that the inversion relation method [18H20] require analyticity assumptions, which cannot be rigorously 
justified, it would be desirable to reproduce the result (|2.42|l by other methods, e.g., by a kind of Bethe Ansatz. 
*See Appendix A for the corresponding expressions in the original variables 



10 



and 

£(^101,02) = -|^(C(0i)+C(02)) 

I </>i-</>2 ,, , , ,N 4>i+(t'2 , , ,N "4 (3.51) 



■)2(^(z + e)|T') i »tfs(J(-' + e)|r') 

Note the leading terms in ()3.48p are independent of the integers nj, entering ()3.46p . Exphcit expres- 
sions for WQ{(j)i, (j)2,ni,n2) and S{4),n) are given in the next Section. 

3.2 Expansion of the star-triangle relation 

In deriving ()3.48p we assumed the normahzation (j2.42p for which the factor IZ in p.Sp is equal to 
one. Substituting ()3.48p and (I3.49|) into the star-triangle relation ()1.5p . one obtains 



(3.52) 



= (r'/r) e--^A(</>)/e We, (^2,^3) VF.-^,,-e3 (m, ng) VF^a (^1, ^2) (l + 0(e)) , 

where 0= ((/<o, </>!, </'2, (As), and 

A^icP) = C{TT-ei\<Pu(Po)+Ci9i+e3\^2,<Po)+CiTT-e3\^3,<Po)+C{<l)o) , (3.53) 

^a('A) = C{ei\<P2,<P3) + C{7^-0i-e3\(l)3,(l>i) + C{9s\(Pi,(l)2) . (3.54) 

Here we used the abbreviated notations S'((/>o, no) = S'(no), Wg{ni, rij) = We{(j)i, (j)j,ni,nj), assuming 
an implicit dependence on the variables (pi. 

Evaluating the integral (|3.52p by the saddle point method one immediately obtains two non- 
trivial identities valid for arbitrary values of , 02 , 03 • The first of these relations (in the leading 
order in e) reads 

•^'k{(p'o^\4>i,4>2,4>3) = -4a (01, 02, 03), (3.55) 

where 0q'^''' is the stationary point of the integral in (I3.52p . i.e., the value of 0o, which solves the 
equation 

~fU — ^ ^(-0 = ° ■ ^^-^^^ 

In the following we will omit the superfix "(cl)" and always assume that 0o = 00*^'^ • To write (j3.56p 
explicitly, define the function 

Using the expressions for £(0 | 0i,02) given ()3.5ip one can write the equation ()3.56p in the product 
form 

^^-0, (00, 0l) ^91+^3 (00, 02) ^^-03 {^0, 03) = 1 • (3.58) 
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(3.59) 



The second relation, which fohows from (|3.52|) (in the order 0{e^)) reads 

= RWe^{cl)2,(p3,n2,n3)WTr-ei-e3{4>i,4'3,ni,n3)Wd.^{(pi,(j)2,ni,n2) , 
where 0o 02 , 03 satisfy the equations (|3.58p and 

i?=(r7r)(^';(0o)/27r)^ . (3.60) 
The quantity ^'^(^o) denote the second derivative of the action ()3.53p at the stationary point (j3.56p . 

3.3 Energy functional 

Using the asymptotics ()3.48p . (I3.49P in ()2.40p and calculating the integral by the saddle point method 
one obtains 

log Z = --£:(0) +0(1) + 0(e). (3.61) 



The energy functional reads 



SicP) = Y^CiOi, I 0i,0,) + ^C(0„) , (3.62) 

{ij) m 



where 



eij = \aijN/{l + NT), (3.63) 
and the variables = {0i, 02, ■ • • , 0m} solve the variational equations 

n^',^^,(0i,0,) = l, i = l,2,...,M (3.64) 

where for any internal site i the product is taken over all edges {ij) meeting at i (the index j 
numerates these edges). The function ^ is defined in (|3.57p . Note a useful sum rule [10] which is a 
corollary of the definition (j2.4ip for lattices of the type described in Section 2.1, 

^ = 27r . (3.65) 

j 

Note that by this property the C(0j) terms, that arise from the first sum in ()3.62p (see p.5ip ). 
exactly cancel out the second sum in p.62p . 

The equation (|3.55p is the classical star-triangle relation. It was introduced in [10] and used there 
to prove the invariance of energy (or action) functionals of the type (|3.62p under deformations of the 
lattice, discussed in Section 2.1, namely under the star-triangular moves (this is the classical analog 
of Baxter's invariance) . Recently the equation (j3.55p was discussed in |30] in connection with 
equation and other integrable equations from [23] (see also [SI])- As noted in [30] the stationarity 
equation (j3.58p could be identified with the so-called three-leg form of the (54-equations, which could 
be written in other equivalent forms. 

Differentiating ()3.55p with respect to 0j, and taking into account (|3.56p . 



d 



(^*(0o, 01, 02, 03) - ^a(0i, 02, 03)) = 0, 3 = 1, 2, 3, (3.66) 
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one obtains another three equivalent forms of (|3.56|) . These equations can be written similarly to 
(I338D . 



^91+03 ('^2, 0o) = ^ei(02,</'3)^e3(02,</'i) , (3.67) 

^7r-6»3(</'3,</'0) = ^ei(03,<A2)^7r-0i-e3(<^3,</'l) , 

All these relations can be brought to the canonical form of the Q4-equation from the Adler-Bobenko- 
Suris list [24J by the substitution 

Then, each of Eqs. (j3.58p and (|3.67p can be reduced to 

r(6'i) {uqUi - U2U3) + T{6i + 6*3) {uoU2 - U1U3) + T{e3) (uqUs - U1U2) 



+T{e{) T{6i + 63) T{93) {U0UIU2US - 1) = . 



(3.69) 



This is an "affine-linear" constraint on four variables uq,ui,U2,U3, which can be solved for each of 
them via simple rational functions of the other three. 



4 Discrete spin solutions of the star-triangle equation 

Before proceeding to the details of the solution of the star-triangle equation (|3.59p . let us explain 
some conceptual modification to the integr ability scheme in this case. We need to consider an 
arbitrary planar graph of the type discussed in Section 2.1 and assign two types of spin variables 
to the lattice sites: classical variables (pi, satisfying equations (|3.64p . and discrete spin variables 
Hi G Zat. Then one needs to solve all equations (|3.64p and determine all 4>i. After that the discrete 
spin model is defined by the weights Wg{(j)i,(l)j,ni,nj) and 5((/)o,no) given by (I4.74p . Previously, 
related ideas were developed in [32t[33]. 

It is important to note that the Baxter's Z-invariance property manifests itself in each order of the 
lower temperature expansion. Namely, the classical star-triangle relation ()3.55p ensures that classical 
action (|3.62p is invariant under the star-triangle moves of the lattice. Likewise the "quantum" star- 
triangle relation (I3.59P and the variational equations (I3.64p ensure this invariance for the partition 
function of the discrete spin model defined by ()4.74p . 



4.1 General solution 



Here we present explicit expressions for the weights solving the star-triangle relation (I3.59P that arise 
in the low-temperature limit, described in the previous section. The value of the factor R given in 
()3.60p corresponds to the normalization of weights inherited from ()2.28p . ()2.42p via the expansion 
(j3.48p . This normalization it not particularly natural. Here we will use the much more conventional 
normalization 

Wei^i,(l)j,0,0) = 1, y(|)^,(|)J ■ (4.70) 
The corresponding value of the factor R is given in (I4.77P below. Define two functions 



rf){(j);n) 



- 0) I r') 



Mh 



n/N 



(4.71) 
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and 



tQ{(t)]n) 



^3(1(0-^) I r') 



n/N „ 



n 



^4(^(fc-i) + 2^(.^-e)|^) 



(4.72) 



They possess the following symmetries 

h;n + N) = r0{-<j), -n) = r0{(j); n), t0{(l); n + N) = te{-(l), -n) 



i), (4.73) 

and r£)((/);0) = tg{4);0) = 1. These functions generalize those used in [7] in connection with the 
Kashiwara-Miwa model. 

Then the weights satisfying (|3.59|) are given by 



We{4>i,(f)j;ni,nj) = rg{(l)i - <f>j;ni - nj)tg{(l)i + <f>j;ni + nj) , 



1 Mw^o + jr<t>o\w) 



(4.74) 

^ M<Po I t') 
Note that the weights are chiral, for generic (pi , 

Wg{(j)i,(j)j;ni,nj) = Wg{(j)j,(j)i;nj,ni), (4.75) 

however 

We{<j)i,4ij;ni,nj) / We{4>i,(j)j;nj,ni) . (4.76) 

The factor R is given by 

R = FeiFesI Fei+e-i , (4.77) 

where all Fg^, Fg^ and -Fgi+^a also implicitly depend on values of (po, 4)2, (f's^ satisfying Eq. (|3.58p . 
The corresponding expressions are given by 



where 

Next 

Pe{(l)i,4>j] 
and finally 
Qeifpo 



= n 



ra=l 



— dii—n H ) 

2 ^^N n' 



n/N 



, (4.78) 



(4.79) 



N-l 

n 

n=0 



MU^ + l) + 2f('^^ - + ^)) MM^ + l) + 2f('^^ + + ^)) 
^1 (f + l) + 2^('/'. - </'j - ^)) ^4(f (n + i) + 2^(<Ai + - e)) 



N-l 

n 

n=l 



TT 



1 



^1 U?™ + - <^0 + 0fc) h^n + — (0 + 



2iV 



TT 



1 



2Af 



vr 



1 



N 2N 



TT 



1 



N 2N 



-1 fc/AT 



(4.80) 



(4.81) 



where k) is an arbitrary permutation of (1, 2, 3). The above expressions involve theta-functions 
with the imaginary period t'/N, i.e., ■i?j(z) = t?j(z| ^), j = 1,2,3,4, and i?'^ = i?2^?3'*?4 denote the 
corresponding theta-constants. Note that 



,(j)j) = Pg{(j)j,<j)i), Qg{(l)i,4ij) = Qg{(j)j,(pi) . 



(4.82) 
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The above solution of the star-triangle equation (|3.59|) contains six continuous parameters: the 
two original spectral parameters ^1,^3, the elliptic period r', Imr' > 0, and the three parameters 
(j)i,(j)2-,<p2,- There is also an integer parameter N > 2, which denotes the number of the discrete spin 
states. 



4.2 Kashiwara-Miwa model 

The simplest solution of (j3.58p is 

(^j=Tr{C + iy), CgZ, ^^£{0,^}, j = 0,l,2,3, (4.83) 

when each factor in p.58p equals to one, ^e{(j)i,(j)j) = 1. This case leads to the Kashiwara-Miwa 
(KM) model [H[71l34j . Since the parameters C and v are the same for all Boltzmann weights, there is 
no need to explicitly display the (/)-dependence in the star-triangle relation (13.590 . It can be written 

as 

E 5(™)(no) H^i™;(m,no) wi'^f,^{n,,n,) W^i™J(n3,no) 

noSZiv (4.84) 

= /?(™)<^')(n2,n3)W^i™;_,^(ni,n3)<™)(ni,n2) . 

where the superscript "KM" stands for the Kashiwara-Miwa model. Introduce functions 

^(KM)/ N _ TT ~\) ~ 2^*^! ^) 7^(km)/ N _ TT ^4(^(fc - \ + v) - ^6 \ jj) 

The expressions ()4.74p give 
wt^\n.,n,)=r^r\n.-n,)t^r\n. + n,+C) , S^^^^n) = ^ (2n + C + .) | ^) ^ ^^^^^^ 

Vjv W4(7rZ/|T'j 

where have we changed the normalization 

W^""''^ (0, 0) = (0/4™^ (0) , (4.87) 

which is more convenient in this case than (j4.70p . The factor R^^'^^ in (|4.84p is obtained by special- 
izing the expression (j4.77p . where Fg is defined in (j4.78p . Moreover one needs to take into account 
the change in normalization (I4.87P in comparison with (I4.70p . In this way one obtains 

p(km) _ p{KM) p(KM) ,p(KM) ( A QQ ) 

^ - ^ei /^ei+ea ' (4.88) 

where 

I JV I I N-2v I 



k=l 



and the symbol [x\ denotes the integer part of x. 

The resulting solution of the star-triangle relation contains three continuous parameters, 61,62 
and t' and three discrete parameters N,C, and = 0, ^. The original Kashiwara-Miwa paper [l] 
only deals with the case = 0. The case v = ^ was considered in [34J- The expression (j4.89p for 
z/ = was obtained in [7] . 
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4.3 Chiral Potts model 

Consider next the limit r' — t- ioo (trigonometric limit). In this limit the Jacobi theta- functions 
i93,i?4 of a real argument x become identically equal to one, i)3{x\\oo) = '(?4(x|ioo) = 1. Therefore all 
functions ^ in (|3.57p depend only on the difference (pi — (pj and thus satisfy the inversion relation: 

cos hichi — 4>j + 6) 1 , ^ 

*^w>*''^^- = — r r m = -^TTTTir^ ■ ^-90 

cos ^[(pi - - e) '^e{(Pj,(Pi) 
Eq. (j3.58p . which determines the value of (po, can now be solved explicitly 

V e+"?^i sin 01 + e+"?^2 sin 6*2 + e+"^3 sin 6*3 y ' 

where 62 = tt — 61 — O^. 

The star-triangle relation (|3.59|) in this case contains only four independent continuous parame- 
ters: two differences (pi — (p2 and (pi — (ps and two spectral parameters 61 and ^3. We will show that 
this star-triangle relation is equivalent to that of the chiral Potts model upon a change of variables. 
The latter has the form (j2.18p containing three rapidity variables p, q and r. Moreover, the Boltz- 
mann weights of the chiral Potts model depend on one additional parameter — the modulus of the 
rapidity curve (see below). 

Each rapidity variable p is represented by a three-vector {xp,yp, fip) specifying a point on the 
algebraic curve, defined by any two of the following equations (the third equation follows from the 
other two) 



x: 



+ <=A:(l + x;^<), kx'; =l-k'^,p'\ ky'; = l-k'fi'; . (4.92) 

where /c is a constant (the modulus of the curve) and k'^ + k''^ = 1. The rapidities q and r are defined 
similarly. The Boltzmann weights of the chiral Pottis model are given by [1] 

(\ n " A; "A; 

where we have assumed the normalization 

Wp,(0) = Wpq{Q) = 1. (4.94) 
These weights satisfy the star-triangle relations (|2.18p where the factor Rpqr is given by 

n(CP) ^ fqr^^fpT^ r{cp) ^ TJ f A^g(l " ^^) (^p " ^^^g) (^^g " ^^j/p) \ ^ 

fj,^^^ ' l}^Wp{xp-u;%){yp-u;%){xp-u:%)) ' ^^""'^ 

Note that for the chiral Potts model the two star-triangle relations in (j2.18p are corollaries of each 
other. 

Remind, that the variables 61,63 and (po,(pi,(p2,(p3, appearing in (|3.59p are constrained by the 
relation (I4.9ip . Moreover the (/>- variables enter only through the differences (pi — (pj. This leaves 
only four independent continuous degrees of freedom. Let us parameterize them via three rapidities 
p, q, r and the modulus k of the curve (I4.92p . 



\Xc{yq)^ [Xpyp)2 



(4.96) 
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One can easily check that with this parametrization the constraint (|4.91|) is satisfied by virtue of the 
equations defining the algebraic curve (|4.92p . With these substitutions the weights (|4.74p exactly 
transform into the weights of the chiral Potts model, given by (j4.93p . For instance, the following 
equalities 

We^{(j)2,(p3,'n2,n3) = W^^^\n2 - ns) , W^-eA4'o,(pi,nQ,ni) = Wg^/\no - ni) , (4.97) 
immediately follow from the relations 

(4.98) 

Vq 

which are simple corolaries of ()4.96p . Note that the normalization (j4.70p precicely corresponds to 
(j4.94p . Similarly one obtains, 

Wg^j^g^{(l)2,<pO,n2,no) = W^r {n2 - TIq) , Wtt-Si -63 ("/"S ,</>!, ^^3 , "l) = (723 - m) , 

W0^{(j)2,(p3,n2,ns) = W^r {n2-ns), WT,-e^{(pi,(po,ni,no) = Wg^. (ni - no) . 

(4.99) 

The Boltzmann weight S^uq) of the central spin in ()3.59p becomes constant, S^uq) = \/^/N , while 
the expressions (j4.78p transform to 



so that the factor (j4.77p exactly reduces to R^qr^ /VN. Thus, there is a precise coincidence with 
the standard form of weights and star-triangle relation for the chiral Potts model [3l[7]. Note also, 
that all square roots appearing in the parametrization (|4.96p . completely cancel out in the final 
expressions for the Boltzmann weights. 

To conclude this Section, let us express the rapidities p, q, r, and modulus k of the curve (j4.92p 
in terms of the O's and (/)'s. It is convenient to define new variables Ai, A2, A3, such that 

01 = A3 - A2, 03 = A2 - Ai . (4.101) 

Eq. (l4.9ip can now be re- written as 

gi^o = ei{0i+02+<^3) C/± = e^'<^isin(A2-A3)+e±'<^2sin(Ai-A3)+e±''^^sin(Ai-A2). (4.102) 

Introduce an additional (dependent) variable Aq by a formula obtained from the previous one by 
interchanging all (j)^s and A's, 

giAo ^ gi(Ai+A2+A3) y± =e±'^isin(</.2-03)+e^'^2sin(</.i-03)+e^'^3sin((/.i-(/)2). (4.103) 

It convenient to define the quantities 



(4.104) 



(Ao + Ai 


-A2 


-A3)/2iV, 


/i = ( 


(po + 4>i- • 


02-.^3)/2A^, 


(Ao + A2 


-Ai 


-A3)/2iV, 


/2 = ( 


(1)0 + (t>2- ' 


01-</.3)/27V, 


(Ao + A3 


-Ai 


-A2)/2iV, 


/3 = ( 


h + 03 - ' 


0i-,/.2)/2iV. 
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Then, inverting (|4.96|) . one obtains 



w2 e 



Xn 



UJ2 e 



i(«3+/3) 



and 



(4.105) 



(4.106) 



= (F+F_)/(C/+C/_) . 

For real i;^'s and A's the quantities V+ and V- are complex conjugate to each other (the same is true 
for C/+ and U-). In this case the last formula can be written as 



k 



e''^^ sinf 



2) + e'-^2 sm{(j)i 



63) + sin( 



e'<^i sin(A2 - A3) + e'^^ sin(Ai - A3) + e'«^3 sin(Ai - A2) 



(4.107) 



It is worth noting an explicit, but rather unexpected, symmetry of the above expressions upon 
interchanging (p's and A's. It would be interesting to understand reasons behind this phenomenon. 

The results of this Section are closely related to the description of the chiral Potts model previ- 
ously obtained by one of us in [33j. 



5 Conclusion 

We have obtained a new solution (jl.ip of the star-triangle relation (|1.5p . which is expressed in terms 
of the elliptic F-function. The solution defines a perfectly physical two-dimensional solvable model 
of statistical mechanics with positive Boltzmann weights. Its partition function is defined in (|2.40|) . 
This is an Ising-type model with continuous real spins < Xj < 27r, which can be regarded as angle 
variables on a circle. The model contains two temperature-like variables p and q; see ()1.6p for a 
definition of the main physical regimes. 

We would like to stress that the formula (11. ip provides a master solution of the star triangle- 
relation in the sense that it contains as special cases all continuous and discrete spin solution of 
this relation that were previously knowiO. In particular, in Sections 4.2 and 4.3 we explicitly 
demonstrate how to obtain the Kashiwara-Miwa and chiral Potts models. We show that these 
models are particular cases of a more general "hybrid" model, which couples a classical integrable 
system involving continuous lattice spins to an Ising-type model of statistical mechanics with discrete 
spins, taking N >2 different values. This hybrid model arises from the general case when one of the 
temperature-like parameters approaches a root of unity, = 1. The required limiting procedure is 
considered in Section 3. To formulate this hybrid model on an arbitrary planar graph (of the type 
discussed in Section 2.1) one needs to assign two types of spin variables to lattice sites: classical 
variables (pi, satisfying equations ()3.64p . and discrete spin variables rii G Ztv- Then one needs to solve 
the equations (I3.64p with fixed boundary conditions and determine the variables <j)i for all internal 
sites of the lattice. After that the discrete spin model is defined by the weights Wgi^cpi, (j)j,ni, rij) and 
S{(po,no) given by (I4.74p . In general, it is a spatially-inhomogeneous model where the Boltzmann 
weight vary across the lattice, since they depend on the variables (pi. Nevertheless, it is an integrable 
model. In particular, its partition function also possesses the Baxter's Z-invariance property of 
Section 2.1 by virtue of the star-triangle relation (I3.59P and the variational equations (I3.64p for the 
classical variables (pi. 

It appears that our master solution (II. ip is also deeply connected with many other important 
models of statistical mechanics, whose consideration goes beyond the scope of this paper. We hope 
to unravel these connections in the future. 

^See footnote 1 on page 2. 
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Appendix A. Energy functional in original variables 

In the original variables ^j,a,T the asymptotic of ()2.28p and ()2.29p is 



logWa(xi,X2) 

log S{xi) 



--/:(a|6,6) + o(i) , 



(A.l) 
(A.2) 



where C{a \ ^1,^2) = | </>i, (^2) and C(^i) = C^cpi), ^, a are related with (p, 9 by ()3.47p and C and C 
have appeared in ()3.48l3.49p . Expressions for C and S are then given by 



and therefore 



^3 f z-jg iVr , 

dz log — j-^j + 

^^3(|(^ + ia)|iVr) 

■k/N 



dz log 



^3{f{z-\a)\NT) 
^3{f{z + \a)\NT) 



(A.3) 



d,iogMi(!^i±i^ + 

^^i(f(ia-^)|A^r) + 



dz log 



n/N 



i?i(f (ia + z) I iVr) 
^9i(f (iQ-z) I iVr) 



vr vr 
Finally, the function C{^i) has the period — and on the interval < < — it is given by 



TT 



(A.4) 



(A.5) 



Note that the quadratic in ^ terms cancel out from the energy functional and, therefore, do not 
contribute to the variational equations (ll.lSp for square lattice. For p = 1 and T = e, the energy £ 
in (jl.lOp is given by 



(ij) {kl) m 



(A.6) 



where the first sum is taken over all horizontal edges {ij), the second over all vertical edges (kl) and 
the third over all internal sites m of the lattice. (11.81). 
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